Evaluation of mixture toxicity literature and chemical space: a data centric approach.

Using the Abstract Sifter tool to analyse PubMed, we reveal published mixture related research most commonly relates to water pollutants, pesticides, environmental pollutants, insecticides, soil pollutants, and chemicals described as persistent, bioaccumulative, and toxic. Furthermore, we discern individual chemicals that also identify as priority chemicals in biomonitoring initiatives and using an ontology-based chemical classification, at the level of the chemical subclass, found these priority chemicals overlap with just 9% of the REACH chemical space.

Quantifying Analogue Suitability for SAR-Based Read-Across Toxicological Assessment.

Structure activity relationship (SAR)-based read-across often is an integral part of toxicological safety assessment, and justification of the prediction presents the most challenging aspect of the approach. It has been established that structural consideration alone is inadequate for selecting analogues and justifying their use, and biological relevance must be incorporated. Here we introduce an approach for considering biological and toxicological related features quantitatively to compute a similarity score that is concordant with suitability for a read-across prediction for systemic toxicity. Fingerprint keys for comparing metabolism, reactivity, and physical chemical properties are presented and used to compare these attributes for 14 case study chemicals each with a list of potential analogues. Within each case study, the sum of these nonstructural similarity scores is consistent with suitability for read-across established using an approach based on expert judgment. Machine learning is applied to determine the contributions from each of the similarity attributes revealing their importance for each structure class. This approach is used to quantify and communicate the differences between a target and a potential analogue as well as rank analogue quality when more than one is relevant. A numerical score with easily interpreted fingerprints increases transparency and consistency among experts, facilitates implementation by others, and ultimately increases chances for regulatory acceptance.

A New Approach Methodology (NAM) Based Assessment of Butylated hydroxytoluene (BHT) for Endocrine Disruption Potential.

Butylated hydroxytoluene (BHT) is a synthetic antioxidant widely used in many industrial sectors. BHT is a well-studied compound for which there are many favorable regulatory decisions. However, a recent opinion by the French Agency for Food, Environmental and Occupational Health and Safety (ANSES) hypothesizes a role for BHT in endocrine disruption (ANSES (2021). This opinion is based on observations in mostly rat studies where changes to thyroid physiology are observed. Enzymatic induction of Cytochrome P450-mediated thyroid hormone catabolism has been proposed as a mechanism for these observations, however, a causal relationship has not been proven. Other evidence proposed in the document includes a read across argument to butylated hydroxyanisole (BHA), another Community Rolling Action Plan (CoRAP)-listed substance with endocrine disruption concerns. We tested the hypothesis that BHT is an endocrine disruptor by using a Next Generation Risk Assessment (NGRA) method. Four different cell lines: A549, HCC1428, HepG2, and MCF7 were treated with BHT and a series of BHT analogs at 5 different concentrations, RNA was isolated from cell extracts and run on the L1000 gene array platform. A toxicogenomics-based assessment was performed by comparing BHT's unique genomic signature to a large external database containing signatures of other compounds (including many known endocrine disruptors) to identify if any endocrine disruption-related modes of action (MoAs) are prevalent among BHT and other compounds with similar genomic signatures. In addition, we performed a toxicogenomics-based structure activity relationship (SAR) assessment of BHT and a series of structurally similar analogs to understand if endocrine disruption is a relevant MoA for chemicals that are considered suitable analogs to BHT using the P&G read across framework (Wu et al., 2010). Neither BHT nor any of its analogs connected to compounds that had endocrine activity for estrogens, androgens, thyroid, or steroidogenesis.

Bioactivity profiling of per- and polyfluoroalkyl substances (PFAS) identifies potential toxicity pathways related to molecular structure.

Per- and polyfluoroalkyl substances (PFAS) are a broad class of hundreds of fluorinated chemicals with environmental health concerns due to their widespread presence and persistence in the environment. Several of these chemicals have been comprehensively studied for experimental toxicity, environmental fate and exposure, and human epidemiology; however, most chemicals have limited or no data available. To inform methods for prioritizing these data-poor chemicals for detailed toxicity studies, we evaluated 142 PFAS using an in vitro screening platform consisting of two multiplexed transactivation assays encompassing 81 diverse transcription factor activities and tested in concentration-response format ranging from 137 nM to 300 µM. Results showed activity for various nuclear receptors, including three known PFAS targets--specifically estrogen receptor alpha and peroxisome proliferator receptors alpha and gamma. We also report activity against the retinoid X receptor beta, the key heterodimeric partner of type II, non-steroidal nuclear receptors. Additional activities were found against the pregnane X receptor, nuclear receptor related-1 protein, and nuclear factor erythroid 2-related factor 2, a sensor of oxidative stress. Using orthogonal assay approaches, we confirmed activity of representative PFAS against several of these targets. Finally, we identified key PFAS structural features associated with nuclear receptor activity that can inform future predictive models for use in prioritizing chemicals for risk assessment and in the design of new structures devoid of biological activity.

Expanded high-throughput screening and chemotype-enrichment analysis of the phase II: e1k ToxCast library for human sodium-iodide symporter (NIS) inhibition.

The sodium-iodide symporter (NIS) mediates the uptake of iodide into the thyroid. Inhibition of NIS function by xenobiotics has been demonstrated to suppress circulating thyroid hormones and perturb related physiological functions. Until recently, few environmental chemicals had been screened for NIS inhibition activity. We previously screened over 1000 chemicals from the ToxCast Phase II (ph1v2 and ph2) libraries using an in vitro radioactive iodide uptake (RAIU) with the hNIS-HEK293T cell line to identify NIS inhibitors. Here, we broaden the chemical space by expanding screening to include the ToxCast e1k library (804 unique chemicals) with initial screening for RAIU at 1 × 10-4 M. Then 209 chemicals demonstrating > 20% RAIU inhibition were further tested in multiple-concentration, parallel RAIU and cell viability assays. This identified 55 chemicals as active, noncytotoxic RAIU inhibitors. Further cytotoxicity-adjusted potency scoring (with NaClO4 having a reference score of 200) revealed five chemicals with moderate to strong RAIU inhibition (scored > 100). These data were combined with our previous PhII screening data to produce binary hit-calls for ~ 1800 unique chemicals (PhII + e1k) with and without cytotoxicity filtering. Results were analyzed with a ToxPrint chemotype-enrichment workflow to identify substructural features significantly enriched in the NIS inhibition hit-call space. We assessed the applicability of enriched PhII chemotypes to prospectively predict NIS inhibition in the e1k dataset. Chemotype enrichments derived for the combined ~ 1800 dataset also identified additional enriched features, as well as chemotypes affiliated with cytotoxicity. These enriched chemotypes provide important new information that can support future data interpretation, structure-activity relationship, chemical use, and regulation.

The Tox21 10K Compound Library: Collaborative Chemistry Advancing Toxicology.

Since 2009, the Tox21 project has screened ∼8500 chemicals in more than 70 high-throughput assays, generating upward of 100 million data points, with all data publicly available through partner websites at the United States Environmental Protection Agency (EPA), National Center for Advancing Translational Sciences (NCATS), and National Toxicology Program (NTP). Underpinning this public effort is the largest compound library ever constructed specifically for improving understanding of the chemical basis of toxicity across research and regulatory domains. Each Tox21 federal partner brought specialized resources and capabilities to the partnership, including three approximately equal-sized compound libraries. All Tox21 data generated to date have resulted from a confluence of ideas, technologies, and expertise used to design, screen, and analyze the Tox21 10K library. The different programmatic objectives of the partners led to three distinct, overlapping compound libraries that, when combined, not only covered a diversity of chemical structures, use-categories, and properties but also incorporated many types of compound replicates. The history of development of the Tox21 "10K" chemical library and data workflows implemented to ensure quality chemical annotations and allow for various reproducibility assessments are described. Cheminformatics profiling demonstrates how the three partner libraries complement one another to expand the reach of each individual library, as reflected in coverage of regulatory lists, predicted toxicity end points, and physicochemical properties. ToxPrint chemotypes (CTs) and enrichment approaches further demonstrate how the combined partner libraries amplify structure-activity patterns that would otherwise not be detected. Finally, CT enrichments are used to probe global patterns of activity in combined ToxCast and Tox21 activity data sets relative to test-set size and chemical versus biological end point diversity, illustrating the power of CT approaches to discern patterns in chemical-activity data sets. These results support a central premise of the Tox21 program: A collaborative merging of programmatically distinct compound libraries would yield greater rewards than could be achieved separately.

Mapping allosteric linkage to channel gating by extracellular domains in the human epithelial sodium channel.

The epithelial sodium channel (ENaC) mediates sodium absorption in lung, kidney, and colon epithelia. Channels in the ENaC/degenerin family possess an extracellular region that senses physicochemical changes in the extracellular milieu and allosterically regulates the channel opening. Proteolytic cleavage activates the ENaC opening, by the removal of specific segments in the finger domains of the α- and γ ENaC-subunits. Cleavage causes perturbations in the extracellular region that propagate to the channel gate. However, it is not known how the channel structure mediates the propagation of activation signals through the extracellular sensing domains. Here, to identify the structure-function determinants that mediate allosteric ENaC activation, we performed MD simulations, thiol modification of residues substituted by cysteine, and voltage-clamp electrophysiology recordings. Our simulations of an ENaC heterotetramer, α1ßα2γ, in the proteolytically cleaved and uncleaved states revealed structural pathways in the α-subunit that are responsible for ENaC proteolytic activation. To validate these findings, we performed site-directed mutagenesis to introduce cysteine substitutions in the extracellular domains of the α-, ß-, and γ ENaC-subunits. Insertion of a cysteine at the α-subunit Glu557 site, predicted to stabilize a closed state of ENaC, inhibited ENaC basal activity and retarded the kinetics of proteolytic activation by 2-fold. Our results suggest that the lower palm domain of αENaC is essential for ENaC activation. In conclusion, our integrated computational and experimental approach suggests key structure-function determinants for ENaC proteolytic activation and points toward a mechanistic model for the allosteric communication in the extracellular domains of the ENaC/degenerin family channels.

Gain-of-Function Mutation W493R in the Epithelial Sodium Channel Allosterically Reconfigures Intersubunit Coupling.

Sodium absorption in epithelial cells is rate-limited by the epithelial sodium channel (ENaC) activity in lung, kidney, and the distal colon. Pathophysiological conditions, such as cystic fibrosis and Liddle syndrome, result from water-electrolyte imbalance partly due to malfunction of ENaC regulation. Because the quaternary structure of ENaC is yet undetermined, the bases of pathologically linked mutations in ENaC subunits α, ß, and γ are largely unknown. Here, we present a structural model of heterotetrameric ENaC α1ßα2γ that is consistent with previous cross-linking results and site-directed mutagenesis experiments. By using this model, we show that the disease-causing mutation αW493R rewires structural dynamics of the intersubunit interfaces α1ß and α2γ. Changes in dynamics can allosterically propagate to the channel gate. We demonstrate that cleavage of the γ-subunit, which is critical for full channel activation, does not mediate activation of ENaC by αW493R. Our molecular dynamics simulations led us to identify a channel-activating electrostatic interaction between α2Arg-493 and γGlu-348 at the α2γ interface. By neutralizing a sodium-binding acidic patch at the α1ß interface, we reduced ENaC activation of αW493R by more than 2-fold. By combining homology modeling, molecular dynamics, cysteine cross-linking, and voltage clamp experiments, we propose a dynamics-driven model for the gain-of-function in ENaC by αW493R. Our integrated computational and experimental approach advances our understanding of structure, dynamics, and function of ENaC in its disease-causing state.